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Abstract 

The ^4 X U{\) flavor model of He, Keum, and Volkas is extended to provide a min- 
imal modification to tribimaximal mixing that accommodates a nonzero reactor angle 
013 ~ 0.1. The sequestering problem is circumvented by forbidding superheavy scales 
and large coupling constants which would otherwise generate sizable RG flows. The 
model is compatible with (but does not require) a stable or metastable dark matter can- 
didate in the form of a complex scalar field with unit charge under a discrete subgroup 
Z4 of the U{1) flavor symmetry. 



I Introduction 



The replication of fermion generations along with small quark mixing angles and large neu- 
trino mixing angles may suggest that the matter content of the Standard Model (SM) trans- 
forms nontrivially under a horizontal flavor symmetry group, Gp- If this flavor group Gp 
were a continuous global symmetry, then spontaneously breaking it would induce massless 
Goldstone bosons in the low-energy spectrum, which are difficult to reconcile with experi- 
ment. One might therefore imagine that either is actually gauged at high energy [I] or 
that the global family symmetry is actually a discrete group [2j rather than a continuous Lie 
group. 

The flavor group 

Gf = A^x U{1) (I.l) 



was originally proposed for lepton masses by E. Ma and G. Rajasekaran [3] and subsequently 
studied by many authors. The U{1) is essentially lepton numberQbut where heavy gauge- 

"'^To avoid generating a Goldstone mode upon spontaneously breaking the U(l), it is convenient to intro- 
duce a new SM-singlet scalar field charged under only a discrete subgroup Z„ with a self-interaction that 
explicitly breaks the U{1). We address this issue in Section VII 
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singlet antineutrinos are assigned U{1) charge zero instead of — 1. 



In this paper we will extend a particular model by He, Keum, and Volkas (HKV) \^ which 
uses the same symmetry group Eq. (I.l). This model proposes that, to leading order, the 



CKM matrix is the identity, and the charged lepton and neutrino mass matrices are "form 
diagonalized" by unitary matrices that result in a tribimaximal neutrino mixing matrix. 



PMNS 



1 1 



TB 



1 



\ 

V2/ 



(1.2) 



independently of the values of the charged lepton and neutrino masses. The main theoretical 
challenge for this model is that the charged fermion mass matrices are obtained using an A4 
triplet vacuum expectation value (VEV) pointing in the direction (1, 1, 1), while the neutrino 
mass matrix is obtained using an A4 triplet VEV pointing in the direction (0,1,0). The 
former leaves a Z3 subgroup of A4 unbroken, while the latter leaves a Z2 subgroup of A4 
unbroken. The two groups do not commute, and so a generic scalar potential will spoil the 
two vacuum alignments. The authors called this the "sequestering problem" and proposed 
one solution based on low-energy supersymmetry. 



In this paper we address the sequestering problem by embedding the HKV model in a frame- 
work which has no superheavy scales, and therefore does not generate large renormalization 
group corrections to quadratic scalar couplings. Thus in this model the sequestering prob- 
lem is solved by showing that it is not a problem to begin with: if we set various quartic 
couplings in the potential to be negligible at some high energy scale, they will remain small 
at low energy simply because the separation of scales involved is not large. 



II Leading order masses and mixing 

The leading order results for quark mixing and neutrino oscillations are the same as in 
the HKV model, so here we will only review the field transformation properties under 
SU{3) X SU{2) X U{1) X X U{l)x and the resulting expressions for the mass matrices. 

The SM fermion fields q, u, d, i, e and the SM-singlet antineutrino fields transform as: 

g~(3,2,+i;3,0), J ~ (3, 1, +|; 1 + 1' + 1", 0) , m ~ (3, 1, -|; 1 + 1' + 1", 0) 
£~(l,2,-i;3,+l), e~(l,l,+l;l + l' + l",-l), iV ~ (1, 1, 0; 3, 0) . (II.l) 

There are four SU{2) x U{1) Higgs doublets, {$a}a = i (where a = 1,2,3 labels the triplet 
representation of A4) and cj): 

$ ~ (1, 2, -i; 3, 0) , ~ (1, 2, -i; 1, +1) . (II.2) 
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If we assume the ^4 VEV alignment ($5) = ("^2) — {^3) = mass matrices for 
the down quarks, up quarks, and charged leptons are of a common form: 

/ yf y'f y'f 

Mf = j^Yjve''-, Yj=\yf uy'^ uj*y'} ] {11.3) 

\yf ^^*y'f (^y'f, 

Here the index / = d,u,e labels the three sectors of charged fermions {Q = +|,— |,— 1, 
respectively) in the SM. 

Thus the mass terms 

Cf = -(/i,/2,/3)M^ +h.c. (II.4) 
are diagonalized by the unitary transformation / = W//mass, where 

e'^f uj*e'^f ue''f\ . (II.5) 



The phase angles a, b, c are chosen to make the physical fermion masses real and positive. 



In Eq. (II. 4) the subscripts 1,2,3 label the triplet representation, while the subscripts 



1, 1", 1' label the different A4 singlet representations [see Eq. (II. 1) . 

In this basis the Yukawa matrix is 

{Ul)^Yf = dmg{V3yfe'''f,V3y'fe'''f,V3y';e"'f) (II.6) 
so the mass of each fermion is an arbitrary coupling constant times a VEV v, just as in the SM. 

The quark mixing matrix is = {ICD^UI. Putting this in the standard form Vq = K-qVcKuPq, 
we find Vckm = I- Thus this model explains why quark mixing is small; the realistic CKM 
angles are presumed to arise from higher order corrections, or from an extension to the the- 
ory. In this work we focus on neutrino mixing and do not pursue quark mixing any further. 
Consequently, from now on we suppress the SU (3) quantum number when expressing group 
transformation properties of the fields. 

As emphasized by HKV, the fermion masses are totally arbitrary while the left-handed mix- 
ing matrices are, up to unphysical phases, fixed purely by group theory, in contrast to the 
alternative possibility that the mixing angles might be related to ratios of masses. 

Ill Neutrino masses 

The leading contribution to neutrino masses arises through the seesaw mechanism. An A4 
triplet of SM-singlet real scalars 

5 ~ (1,1,0; 3,0) (III.l) 
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with an A4 vacuum alignment (Si) = {S3) = 0, (5*2) = f 5 7^ (with vs having either sign) 
gives a Majorana mass matrix for N: 



M. 



N 




X 



VnVs 



(III.2) 



Here m^v is a bare Majorana mass term, and i/n is the Yukawa couphng of S to {NN)^^. By 
rephasing the N fields we can take m^v real and positive without loss of generality. In this 
case, the parameter x is in general a complex number with an undetermined phase The 
Dirac mass term arises from -Cyuk ~ yu4>'^^N and is proportional to the identity: 



Md = ttidI , mo 



(III.3) 



where :^'?^</) e*^'* = (0*^). In the seesaw limit the light neutrino Majorana mass matrix is 

1 



m 



D 



rriN 1 — 




(III.4) 



This is diagonalized in the form (W^) M^W^ = = diag(mi, m2, ^3) by the unitary matrix 

_en/2N 




(111.5) 



where the phase angles a, /5, 7 are chosen such that the physical neutrino masses rrii are real 
and positive: 

(III.6) 



mi 



\l + x\ ' 



1712 = , 



\1 — x\ 



where m^^ = |m|,/mAr|. Neutrino oscillations constrain the mass-squared splittings Am^tm ~ 
m| — 1712 ^"^^ ^"^soi ~ "^2 ~ but not the overall scale rriiy. Hence to compare with 



experiment one computes the ratio 



A7^_ r 5.05 - 6.40 iAmi^>0) 
^ = ^' ' " 14.93 - 6.29 (Aml^ < 0) ^^"'^ ^ = 



(III.7) 



ordering mi < m2 < {Aml^^ > 0) and the "inverted" mass ordering < rrii < m2 



This model has been shown to be compatible with Eq. (HI. 7) for both the "normal" mass 
iri] 

The neutrino mixing matrix is 



1 ei(t-ce-f ) 



J_pi(f-«e) 

1 e^d-be+f) 

V 3 





v2 



1 pj(fe-Ce 



J_pi( 



e^2 



■f) 



(III.8) 
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Thus, putting this in the standard form 14 = JC^Vpum'Pu reproduces the tribimaximal mix- 
ing matrix of Eq. (1.2): I^mns = Vtb- 



The recent measurement of a nonzero ~ 0.1 shows that this I^mns cannot be exact 
at low energy. Just as for the leading order prediction Vckm ~ I, this model predicts that 
^MNS ~ at leading order, subject to small but nontrivial corrections. 

In this work we will use the fact that the Z3 and Z2 subgroups of A4 do not commute 
to generate a realistic neutrino mixing matrix. We propose a model in which the Z2 remains 
a good approximate symmetry of the neutrino sector, but corrections due to the breaking of 
Z3 feed into the neutrino mass matrix and generate the nonzero V^s required to fit data. 



IV Scalar potential 

As explained previously, we would like the charged lepton masses to arise from a VEV 
($°) oc (1, 1, 1), which leaves invariant a Z3 subgroup of A4, and we would like the neutrino 
masses to arise from a VEV (S) oc (0,1,0), which leaves invariant a Z2 subgroup of A4. 
These two subgroups do not commute, and so the most general potential does not admit this 
configuration. This is called the sequestering problem. 



In the absence of a more elaborate framework, such as extra dimensions or supersymmetry, 
it is unlikely that this problem can be solved completely, but it can be sufficiently mitigated 
in a technically natural way. The general philosophy is to declare by fiat that certain cou- 
plings are small at a particular energy scale, and then to construct a model for which these 
couplings do not fiow substantially under the renormalization groupj^ 



To make sure no superheavy scales are required, we need to bring down the scale of m^v from 
the usual seesaw scale; for example, if rriN ~ TeV then the neutrino mass scale rriy ~ 10^^ 
eV can be obtained if rriD ~ rrie ~ 10~^ MeV In contrast, the simplest seesaw model has 
rriN Mqut ~ 10^^ TeV and m^) ~ f ~ 10~^ TeV. Since the Dirac mass mo arises from the 
Higgs field (p [see Eq.(III.3)], we want to make the VEV of small while keeping its mass 
above the weak scale. 



The problem of having heavy Higgs doublets with small VEVs was proposed and solved 
in the context of a different model [7J. In the present case we introduce a Gsm x A4-invariant 
complex scalar field Y which carries —2 units of U{l)x charge: 

y~ (1,0; 1,-2) (IV.l) 

and a second y44-invariant Higgs doublet 

0'~(2,+i;l,+l) (IV.2) 

^This is not an unfamiliar concept. In the SM, we declare by fiat that the electron Yukawa coupling is a 
dimensionless number of order 10~^. The reason we are able to do so consistently is that a small value for 
this coupling is protected by chiral symmetry from receiving large corrections. In other words, we fix a small 
value for a dimensionless coupling and show that it does not flow significantly under RG. 
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with the same U{l)x charge as (p but with opposite hypercharge [see Eq. (II. 2)]. The 
dimension-2 operator (piS^^cp'j ~ (1, 0; 1, +2) is invariant under Gsm x ^4 and can couple 
to Y to form an invariant dimension-3 operator. Consider the following scalar potential: 

ViY, 0, 0') = M^yty + M20t0 + M|,0't0' + Ay (Fty)2 + A^K0'V')' 

-ifiY^ie'^^'j + h.c.) + ... (IV.3) 

Here the ellipses contain all renormalizable terms consistent with the symmetry G^m x Gp 
whose roles are assumed subdominant in what follows. We take My < and M|, < as 
usual, but M| > 0. When Y and 0'° obtain VEVs, they induce a VEV for 0° according to 
the relation 

•V^^.... (.V.4, 

where = \/2{(lP) and v^f,/ = \/2{(j)'^) are assumed real and positive for simplicity. Since 
there is no large separation of scales in this model, it is a self-consistent (though arbitrary) 



assumption to declare that the contributions denoted by in Eq. (IV.4), which arise from 
quartic interactions in the Lagrangian, are subleading. If we can drop those contributions, 
then 1(0°) I ~ Thus the VEV is inversely proportional to the square of the mass. For 

example, if /i ~ v^' ~ 10~^ TeV and (Y) ~ TeV, then 

^— ) X 10^ TeV . (IV.5) 

In this way we will assume that the SU{2) x U{1) doublet is the heaviest particle in the 
theory and thereby specify the values of all coupling constants at the scale M. = ~ 
10^ TeV. This is much lower than the superheavy scales that often appear, e.g. Mqut ~ 10^^ 
TeV, and so the renormalization group running to scales far below Ml will not be significant, 
and the non-commuting Z3 and Z2 subgroups of can remain sequestered up to small 
corrections. 



V Minimal modification to tribimaximal mixing 

Consider the case for which Z2 = {/, T2} is unbroken but Z3 = {/, c, a} is broken by small 
perturbations that feed into the neutrino mass matrix. Since ui and z/3 flip sign under Z2, the 
mass terms z/iz/2 and Z/2Z/3 are still approximately zero. However, now that c : 2^1,2,3 — > ^^2,3,1 
^3,1,2 is broken, the terms z/iz/i and 1/31/3 are no longer forced to be equal. The neutrino mass 
matrix in the G^-basis is now of the form [8j 

/a-e /3 \ 
M, = I 7 . (V.l) 
\ (3 a + ej 

If we assume that, to a good approximation, the charged lepton mass matrix remains di- 
agonalized by the unitary matrix studied previousl}]^ then the neutrino mixing matrix is a 

^Note that this Z2 leaves the A4 singlets 1, 1', 1" invariant. This can be seen from decomposing the 
product of two A4 triplets into irreducible representations: ii u,v ^ 3, then the three singlets are of the form 
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one-parameter deviation from tribimaximal mixing with the middle column oc (1, 1, 1) un- 
changed, which implies that a nonzero reactor angle 613 is generated. To leading order in e, 
we have 

iKsl ~ ■ (V-2) 



To search for a suitable high-energy model that will result in the desired perturbation, we 
first study the low-energy effective interactions which are invariant under Gsm x ^4 between 
$ ~ (2,—^; 3,0) and the lepton doublets i ~ (2, — ^;3,+l) . Consider the dimension-5 
interactions 

^dim-5 = CiCi + CvOy + Cv>Ov> + C3O3 + h.C. (V.3) 

where: 

Ci = + ^\^\ + + ^2^2 + hh) , (V.4) 

Oy = ($1 $1 + uj*^l^\ + uj $^$^) (^1^1 + uj £2^2 + 00* kk) , (V.5) 
Ov = +w$5$5 + w*$^$^)(£i£i + a;*£24 + w£34) , (V.6) 

O3 = {<l>l<l>l<l>l<^l<^m) ■ {i2h,hhJii2) . (V.7) 

These operators are invariant under Gsu x ^4 but break U{l)x by two units. If we couple 
any of these to a Gsm x A4-invariant complex scalar Y with U{l)x charge —2, then these 
operators are promoted to dimension-6, and the coefficients Cj are proportional to (Y). 

Consider the case C3 = and ($5) = ("^3) 7^ ('^2)- Defining for convenience the overall 
scale of the VEVs and the small splitting as 

a = 3($;^)2 , 6 = (<l>°^)2 - (<l>?^)2 (V.8) 

we find: 

ciOi + ci'Oy + ci'/Ci" = [aci + (ci + ijj*ci> + u ci")6] ii£i 

+ [aci + (ci + ci' + ci")5] ^2^ 

+ [aci + (ci + w cr + uj*cr>)5] isi^ . (V.9) 
The perturbation by 6 recovers a neutrino mass matrix in the form of Eq. (V.l) with: 



i^{cv-cv')6. (V.IO) 

Thus we seek a model in which ci ~ ci> ~ ci" 7^ and C3 ~ 0, and in which ($°) oc (1, 1, 1) is 
perturbed in the form oc (x, 0, x). In order to have the flavor-depende nlp| p erturbation 



by S not be overwhelmed by the flavor-independent overall shift by a [Eq. (V.8)], we could 
attempt to construct a model in which ci ^ cr ~ 



uiVi +ujP'^U2V2 +a;^^U3W3 with pi,P2 — 0, 1, 2. Since wiWi, M21'2j and M3W3 do not transform under Z2, neither 
does any hnear combination of them. 

^We are not working in the e,^,T basis; we use the word "flavor" here to mean that the perturbation 
depends on the A4 triplet label 1, 2, 3. This usage is consistent with the modern understanding of "flavor" as 
meaning "copy of representation of Gsm" ~ in other words, copies of identical covariant derivatives - rather 
than having any a priori relation to the charged lepton mass basis. 
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VI Model for nonzero V^^ 



One simple way to realize the effective operators from the previous section is to introduce 
SU{2) X f/(l) triplets, A, A', A", whose transformation properties under SU{2) x f/(l) x 
^4 X U{l)x areQ 



A~(3,-l;l,+2), A' ~ (3,-1;!', +2), A" ~ (3, -1; 1", +2) . 



(VI.2) 



The field At can couple to (££)i, the field A't ~ (1')* = 1" to {U)r, and the field A"t to 
The product {A^Y^{^i^j)i is invariant unde r SU{ 2) x [/(I) x A4, but it carries charge 
+2 under U{l)x- Fortunately the field Y in Eq. (IV.l) has just the required property to 
allow the dimension-4 interactions Y'^A'^{^^)i, FtA't($$)^,, and Y^A'^{^^)i/i. 



First consider the following Lagrangian for A: 

jCa = -Mltr(AtA) + {{A^y^ [Ai£^ij)l + B*Y\^i^j)i] + h.c.} 



(VI.3) 



Here A and B are couplings, and the ellipses denote other terms which are not important for 
the present discussion. We imagine thaf\v < MA<M^r~. 10^ TeV. At scales far below Ma, 
the path integral for A can be performed in the Gaussian approximation, which removes A 
from the low-energy spectrum[^ and leaves behind an effective Lagrangian: 



Mi 



[\A\\i^iM^^'^^')i + \B\^\Y\\<^.<l>M'^^''^^')i + {ABY{<l>^'<l>^^U{idj)i + h.c.)] 

(VI.4) 



The last term in Eq. (VI.4) is an explicit realization of the first term in Eq. (V.3), with the 
coefficient 

(VL5) 



ci =AB 



Mi 



This procedure can be repeated for A' and A" to obtain the second and third operators in 



Eq. (V.3) with coefficients 



cr 



A'B'^ 

Ml, 



Cyi 



A"B" 



(VI.6) 



^We express the SU{2) components of the complex triplet A as a 2-by-2 symmetric matrix: 



A., = 



AO 



A— 



(VI.l) 



where the superscript denotes the electric charge. 

^The notation "a < 6" is to be interpreted as "a can be somewhat less than or possibly comparable to b" . 
Also, one could instead take < Ma, which would amount to pushing the scale at which Gf is imposed 
up to = M^. As long as this is still substantially smaller than the usual AfcuTi the main results of this 
paper remain unchanged. See Appendix [C] for a toy model. 

^Instead of working in the Wilsonian framework of performing the path integral one step at a time, the 
reader may prefe r to s tudy the scalar potential with all fields left in the spectrum. In that case the in teraction 
A'fy$$ of Eq. (VI.3) induces a small VEV for A", which will have the same effect as Eq. (VI. 5). In the 



neutrino literature this is known as the Type-II seesaw mechanism [51 [TO] • 
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To leading order there is no term involving A, A', or A" that couples to the A4 triplet 
(^2^3, ^3^1, ^1^2), so we have C3 = 0. By taking Ma large compared to Ma' and Ma", or by 
taking AB small compared to A'B' and A"B", we can make ci ^ cy, cyi. 



Now that we have generated the Cj, we need to generate the Zs-breaking perturbation 5. 
The perturbation away from ($°) oc (1, 1, 1) can be achieved through a cubic interaction 



Kubic($, S) = ft, 5- ($t$)3^ + f,, 5- ($t$)3^ 



(VI.7) 



If we assume that ($°) = -^Va are real, then the two terms contribute to the classical 
potential through one combined invariant: 



(VI.8) 



where /i = /ii + /i2- This is the case studied in Appendix |Bj We find that vi = V3 V2 is a 



solution, with 



(0) 
A<j) 



(VI.9) 



where v^g^ is the leading order value for (5*2). Using Eq. (V.2), we fincQ 

ki' — ci"\ 



e3 



1/(^)1 



~ (0) 



2v2 X^rrii, 



(VI.IO) 



where f{x) = x/{l — x^) and x = yNvf^ /rriN- The various parameters in Eq. (VI.IO) can be 
arranged to accommodate the measurement V^3 ~ 0.1. 



VII Goldstone bosons and dark matter 

Now that we have achieved a realistic neutrino mixing matrix, let us discuss breaking the 
U{l)x to a discrete subgroup to avoid generating a massless Goldstone boson after flavor 
symmetry breaking. Consider the Gsm x A4-invariant scalar Y, which carries charge —2 under 
U{l)x- Instead of the full U{l)x group, let us instead insist only on the discrete subgroup 
: Y — > e*'^+^)^'^/'^y = — y. Then the most general renormalizable scalar potential for Y is 

Vy = M^Y^Y + (im^y2 ^ ^ \y(Y^Yf + {\\'yY^ + h.c.) . (VII.l) 

The terms Y"^ and Y'^ explicitly break U{l)x and remove the unwanted Goldstone mode from 
the low energy spectrum. The U{l)x representation ip m for the other fields should now 
be understood as the transformation property 

SO that all fields with m = 1 get multiplied by a factor of i. 



The Ci have mass dimension —1, so Ve3 is dimensionless as required. 
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Since various aspects of the model [Eq. ( |IV.4 ) and Eqs. (VI. 5), (VI. 6)] depend crucially 



on a nonzero {Y), we take the parameters in Eq. (VII.l) to be such that KeY or Imy (or 
both) obtains a nonzero VEV. There is no symmetry principle that conserves F-number, so 
this SM-singlet does not provide a good candidate for dark matter. 



Given this, one might be motivated to consider a Gsm x A4-invariant complex scalar X 
which carries unit charge under Z^: 

X~ (1,0;1,+1) . (VII.3) 
The scalar potential involving only X and Y is V{X, Y) = Vx + Vy + Vx^^^ + V^^^^^^, where 

Vx = MlX^X + \x{X^Xf + aX'xX^ + h.c.) , (VII.4) 



Vy was given in Eq. (VII.l), and: 

lA-^bic = [i(py + p'Y^)X^ + h.c] , (VII.5) 
^^uartic ^ XxyX^X Y^Y + {\}^xyX^ X Y^ + h.c) . (VII.6) 

If we take Mx positive so that Vx does not induce a VEV for X, then even when Y picks a 
VEV and breaks Z^, the potential V{X, Y) leaves a discrete X-parity Px '■ X — )• —X unbro- 
ken. Thus, if the (X, y)-sector of the theory were secluded from the rest of the Lagrangian, 
the lighter of Re(X) and Im(X) would be stable and provide a candidate for dark matter 

mm- 

Now consider the interaction of X with the other fields. In addition to the obvious "Higgs- 
portal" interactions such as X^X$^$ which conserve X-number, there are quartic interac- 
tions which are linear in X: 

Vp^ = KiX(f)'^'$i- S + K2X(f)'^'eij$^' ■ S + h.c. (VII.7) 

If we suppose that Ki or K2 is small but nonzero, then Px is broken and X can decay. In 
view of the scales in our model, one might suppose that Mx ^ TeV so that X decays either 
to the three-body final state 0*$S' or (p'^S, or if such a decay is kinematically inaccessible, 
then X can decay to whatever final states can be reached through the propagators of virtual 
scalars, e.g. TeV-scale neutrinos or SM charged leptons. For the sake of a rough calculation 
suppose that Mx 3> M^/ + M$ + Ms and estimate the width (inverse lifetime) of X as 

Tx' ^Tx- j^Mx (VII.8) 

where k is either ki or K2, and (27r)" is a phase space factor for n outgoing particles. Recent 
measurements of cosmic ray spectra [131 [E] may indicate the presence of late-decaying dark 
matter, which would require [TB] a lifetime much larger than the age of the universe: Tmin ~ 
10^''s ~ 10^^ TeV~^, where as usual we work in units h = 1. Thus, if we would like to 
interpret the cosmic ray anomalies in terms of our metastable X particle, we need: 

rmin \Mx J \ K J 
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Thus, as is the usual case with a decay rate resulting from renormalizable interactions, we 
require extremely small coupling constants, k < 10~^^, if we insist on a DM interpretation of 
the data. It is worth recalling [see below Eq. (VII. 6)] that if we were to impose Px '■ X — )■ —X 
on the model, then we would have Ki = K2 = 0. Taking this argument in the other direction, 
we note that in the limit 2 — )■ the transformation Px '■ X — )■ —X emerges as a symmetry 
of the Lagrangian. Thus any value of k in Eq. (VII. 9), no matter how small, is technically 
natural in the sense of 't Hooft. In order to satisfy Eq. (VII. 9) we would need to fix ki and 
K2 to be small dimensionless numbers, but we would not need to engineer any unnatural 
cancellation between large numbers. 



VIII Discussion 

We have embedded the HKV model into a framework whose largest scale Ai is much smaller 
than the traditional seesaw scale: A^/Mgut ~ 10~^°. The main requirements of such a 
scenario are an additional y44-singlet Higgs doublet, 0', and an SM-singlet complex scalar, Y , 
which carry U{l)x charge +1 and —2, respectively [Eqs. (IV.l) and (IV. 2)]. Then the scalar 
potential can be arranged such that the small Dirac mass for the neutrino can be generated 



by a Higgs doublet whose VEV is small but whose mass is large [Eq. (IV.4)]. 



This ameliorates the sequestering problem that typically afflicts models of this type for the 
simple reason that the separation between the scale at which = A4 x U{l)x is imposed 
and the scale of neutrino masses is not sufficiently vast to induce large renormalization group 
corrections. This is an alternative to the supersymmetric and brane-world approaches men- 
tioned by HKV and studied by various authors. 



The recent measurement of Vez ~ 0.1 is accommodated by a minimal modification to 
tribimaximal mixing due to A4-singlet Higgs triplets which communicate the breaking of 
Z3 = {/,c, a} into the neutrino sector while leaving the Z2 = {11^2} essentially unbroken 



(see Section VI). 



The desire to remove an unwanted Goldstone boson due to breaking the continuous global 
U{l)x led us to insist only on the discrete subgroup < t/(l)x, and the Lagrangian for 
the SM-singlet Y explicitly implements the distinction between the two symmetry groups 



[Eq. (VII.l)]. This then motivated (but did not require) the introduction of a second SM- 



singlet complex scalar field X with half the Z4 charge of F, which can be a candidate for 



decaying dark matter [Eqs. (VII. 3) and (VII. 7)]. 



We conclude by mentioning other possibilities for neutrino mixing in this framework; the 
one-parameter deviation from tribimaximal mixing studied here is by no means the only 
option [16]. For example, the study of Appendix |b] might lead us to consider the case where 
Z3 = {/, c, a} remains a good symmetry of the neutrino mass matrix, while the breaking 
of Z2 = {/, ^2} becomes significant. In this case, the neutrino mass matrix in the Gi?-basis 
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becomes 

(a ei2 /3 
ei2 7 ^23 I (VIII. 1) 

P £23 « 

with small parameters Eij. If W£ remains approximately unchanged, then Ves = [(W|)^W|^]e3 7^ 
must arise by modifying the third column of U^. For example, 



M J =(a±/3) +(£i2±£23) 1 . (VIII.2) 






Note that if = +^23? then the third column of remains unchanged from its leading 
order value, which is not what we want. 
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A A4 invariant potential for a real triplet 

The scalar potential for a Higgs field, (p, transforming as an A4 triplet was studied previously 
by one of the authors [1^. Invariance under SU{2) x U{1) implies that this potential is a 
function only of the combination flfb, where a, b label the A4 components, and it has one 
quadratic invariant and five quartic invariants. 



Here we study the scalar potential for a real scalar field S transforming as a triplet un- 
der A4. Since 5* is real, the five quartic invariants for ip collapse to only two independent 
terms. However, since we do not impose a parity that flips the sign of 5* (the real-valued 
analog of a f/ (1) symmetry), there is also a nontrivial cubic invariant from the group theoretic 
property 3 x 3 x 3 = 3 x (H- 1' + 1" + 3^ + 3^) = 3 x 3^ ... = 1 where the 1 denotes 
the invariant S1S2SS. (Note that 3a = since SaSb = SbSa-) 

The most general renormalizable, y44-invariant scalar potential for ~ 3 is 

ViS) = \M' Si + ^^S,S2S, + \\ Si] + ^A' [{S.S^f + (52^3)^ + {S:,S,f] . 

a = l \a = l J 

(A.l) 

We study three|^ potential vacua of interest: E : {{Si) = {S2) = {S3) = S ^ 0}, U : {{Si) = 
S^O, {S2) = {Ss) = 0}, and P : {{Si) = {S2) = S ^ 0, {S3) = 0}. 

^The case (5*2) = {S3) — —{Si) was studied by E. Ma in the context of another model based on A4 
symmetry |18j . 
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In the vacuum E we have = —(3 A + 2X')S'^ — fiS, and: 



v\^ = -U\ + lx')s'-ys' 



and 



dSadSt 



A A 
A A 
A A 



(A.2) 



where = 2A5'^ — fiS , A = 2(A + A')5'^ + fiS. The notation V\^ denotes the value of the 
potential evaluated at the minimum for the vacuum E. Diagonalizing the matrix of second 
derivatives gives V"^ = diag(m^ + 2A, — A, — A). The requirement — A > implies 



In the vacuum U we have = —XS'^, and: 
V\^ = -\XS^ and 



dSadSt 



ml 
A 
A 



(A.3) 



where ttt-q = 2AS'^, = (A' — X)S'^, and A = fiS. The matrix of mass-squared eigenval- 
ues is Vfj = diag(mQ,m^ — 2A,m^ + 2A). The requirement — 2 A > is satisfied if 
liS < |(A' - X)S\ 

When studying the putative vacuum P we find an important difference from the previously 
studied Higgs potential. Minimizing the potential with respect to 5*3 results in the equation 

dV 



fxS' = 0. 



(A.4) 



If /i 7^ then this equation cannot be satisfied. Thus, unlike the case of a complex SU{2) x 
[/(l)-invariant Higgs potential, a real scalar triplet does not admit the vacuum P even as a 
local minimum 



B Interactions between a Higgs field and a real triplet 

Now we study in detail whether the VEVs {S) oc (1,0,0) and {ip) oc (1, 1, 1) can hold to a 
good approximation in a Lagrangian based on x Z2. Let ~ 3_ and S* ~ 3+ be two real 
triplets of A4, where the subscript indicates the parity under Z2. The real field (p which is 
odd under Z2 serves as a toy model for a complex field $ which is charged under SU (2) x U{1). 

The potential isV = Vs + V^ + V^^^'"" + 1^5'7''*'^ where 
Vs = lMlJ23l + f^sSiS2Ss + lXslj2sl] + ^A'^ [(^1^2)' + (^2^3)' + (^3^1)'] (B.l) 

a = l \a = l / 

= IK J2^" + I^AJ2v"] + IK [(^1^2)^ + (^2^3)^ + i^s^i)'] (B.2) 

a = 1 \a = 1 / 

^"in the SU{2) x U{1) case P is a local mininium but never the global one [T7] . 
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describes each sector, and 

^5^'" = KSl^2Vz + S^VZVI + 53<^1<^2) (B.3) 

ysT^'"" = i^o(^Sl)(^ipl) + ni{Slipl+ cyclic) + K2{^^^ 

a b 

(B.4) 

describe the interactions between the two sectors. 

We would hke to find a perturbative solution around the configurations S oc (1,0,0) and 
(p oc (1, 1, 1). The general case can be parametrized as 

^„ = ^f v.. = v^i") + ev^i^) (B.5) 

where e is a small parameter. We will first assume that y^^^''*"^ can be dropped, or in other 
words that the are at most C(e^) and thus negligible. We will assume that V^^*^"^ = eVi 
is itself a perturbation by defining /i = eflo with /io a parameter of 0{fis)- The potential 
Vq = Vs + is to be treated as the leading order potential of 0{e^), with any possible 
renormalization group corrections to its parameters being of O(e^) or higher. 

The goal is to solve the equations g|-(Vo + eVi) = and g^(Vo + eVi) = to 0{e). In 

the particular leading order vacuum of interest, Sj^*^ oc (1,0,0) and ip^'^ oc (1,1,1), this 
procedure gives the result 

Sa = ^0 + eSi'^ , = Q y.0 + e^^^^ (B.6) 

where: 

_ 1 f ~ nil) _ nW _ 1 f "PoV r, m 7^ 



'0 



(1) _ A^ + A'^ f ~ (1) _ (1) _ , -^'^ ^0 



If we take v^o/'S'o ^ 1, then the corrections to Sa oc (1,0,0) can be ignored safely, but those 
to ipa oc (1, 1, 1) are important. This is precisely the situation we seek in the main text: 
the Z3 "charged lepton" subgroup of is broken, but the Z2 "neutrino" subgroup of ^4 
remains conserved to a good approximation. The effect we are talking about is of order 
Si'^/i^? - y^r^) ~ MSof, e.g. ~ (1/4)3 ^ 



A convenient way to parametrize this breaking of Z3 < A4 is: 

_~ 'S'o 
A 



S^Uv'2-Vi)=efio^ + Ois'). (B.9) 



This is the equation we use to obtain Eq. (VI. 9) in the main text. 
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C RG for toy model coupled to fermions 



Consider a simple toy model of two real scalars 0i, 02 coupled to "quarks" q and "antiquarks" 
g, with a color gauge group SU (N) under which q and q transform as and N*, respectively. 
(The (pi are color-neutral, and the gluons are denoted by {6-;^}^=-/.) We impose the reflection 
symmetrjj^ 

Z2 : 01 -01, g -g, 02 ^ +02, g +g (C.l) 
so that 01 has a Yukawa interaction with gg whereas 02 does not: 

>^^Yuk = -y<Piqq + h.c. (C.2) 

The most general renormalizable potential for 0i and 02 consistent with Z2 symmetry is: 

V^(0i, 02) = \MlcPi + \mUI + M12 0?02 + lM\ + hMl + iAi20?0^ . (C.3) 

We are interested in the renormalization group equations (RGE) for the Lagrangian 

Ctoy{G, 01, 02, g, g) = Ci,m{G, 01, 02, g, g) + £Yuk(0i, g, g) - 1^(0i, (C.4) 

in the deep UV regime, far above the scalar masses Mi,M2 and the dimensionful cubic 
coupling M12. (Since we are not interested in low-energy phenomena such as spontaneous 
symmetry breaking, we take and M| both positive.) 

The issue we would like to address is the following. Let m denote a "low-energy" mass 
scale m ~ Mi ~ M2 ~ M12, and let Ml denote the UV cutoff of the field theory. Suppose we 
(arbitrarily) fix the parameters of the theory at the scale Ai to satisfy 

Ai(A^) ~ A2(7W) -0(1) , |Ai2(7W)|<l. (C.5) 

Our goal is then to find the maximum ratio M./m such that the "off-diagonal" coupling A12 re- 
mains parametrically smaller than the "diagonal" couplings Ai and A2 at energy ~ m <^ A^. 

This toy model is a simple subcase of a general class of theories studied by Cheng, Eichten, 
and Li [T9]. Following their notation, we express our Lagrangian a^^ 

£ = -\G%G''^''-\d,cP,d^cP, + tq^^q + tq^pq-{qKqcP, + /i.c.)-K(0) , (C.6) 
'^{4') = ^ fijki4>i(t^j4>k(pe + lower-dimension terms . (C.7) 

Here Gjjj, = d^Gl — dyG'^^+gG"'^'^G^^G1 is the gluon field strength, and ^ = d-^^D^ is the covari- 
ant derivative for each 2-component fermion: D^q = {dfj_ — igG'^^t'^)q and D^g = {df^ — igGp:"')q 

with r = -(r)*. 

^-"^This Z2 serves as a toy model for SU{2) gauge symmetry, under which Higgs fields and left-handed 

quarks transform, but under which SM-singlet scalars ("02") and left-handed antiquarks are neutral. 

-'^^In contrast to the more general case studied in the reference, but similarly to the realistic case in the SM 
and many of its extensions, here our scalars 0i are not charged under the color group. Also here we need 
only one flavor of quark, in which case hi is a single number (for each i) rather than a matrix. 
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The lowest-ordei RGE for the couphngs g,hi, and fijke are (P = IGir'^d/dt, t = In 



for the scale /i and an arbitrary reference point /io): 



Mo 



(C.8) 



— 3(2t hit +5*2/^1)5' + 2887r2 fijkefjkimhm 



(C.9) 



T-^ fijkl fijmnfmnkl ~l~ fikmnfmnjl ~l~ fi£mnfmnjk ~l~ ^^^{hihm) f mjkl 

-2tr{hihj{hk, he} + hihk{hj, hi} + hihe{hj, hk}) 

acd^bcd 3^6"^^ f^f^ 



Here we have defined the group theoretic quantities C C 
'3 



tr(t"t^) = 3^5°''^ related to the color representation of the fermions 



1 

2' 



14 



(C.IO) 



S'f /, and 



For SU{N), we have 



For our case [Eqs. (C.2) and (C.3)] we have: 



hi 



/nil 
/1122 
/1112 



Al , /2222 = 

/1212 = /1221 
/1121 = /1211 



A2 



/2121 
/2111 



/2112 
/2221 



/2211 
/2212 



A12 5 
/2122 



2221 



Then Eq. (C.9) simplifies to 



and Eq. (C.IO) produces three reasonably simple coupled equations: 



V\i = 3{Xl + XI2) + Sy^Ai - 12y^ 



VX2 = 3{Xl + XI2) 



VX 
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(Ai + A2)Ai2 + 4A?2 + 8y2Ai2 



(C.ll) 
(C.12) 
(C.13) 
(C.14) 



(C.15) 



(C.16) 
(C.17) 
(C.18) 



Before solving the system numerically, let us first obtain some intuition for how the "off- 
diagonal" coupling A12 flows. The whole point of our study is to consider the case for which 
A12 is parametrically smaller than the other couplings, so suppose we can drop terms of 



0(A^2)- If 'we drop the 0(A^2) term in Eq. (C.18), then the flow equation for A12 is of the 
form I'Ai2 = J-'(/i)Ai2(yu), with = Ai + A2 + 4Ai2 + 8?/^, and the equation can be integrated: 



Ai2(m) ^ XuiM) e- ^^^"'^ . 



(C.19) 



^^As pointed out in the reference, the last term of Eq. ( |C.9 ) is of "lowest-order" even though it arises from 
a 2-loop diagram. 

^''Our two-component quarks are coupled in a "vectorlike" fashion to the gluons, so it is convenient to 
package them into 4-component Dirac spinors Q — {q,q^), which here transform under the defining A^- 
dimensional representation of SU{N).] 
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Thus in the hmit X^iM.) — ?• 0, the couphng Ai2(/i) gets renormahzed multiphcatively, so 
perhaps contrary to intuition, we should expect that the assumption Ai2(/i) ^ 1 remains 
vahd even for yU ~ m ^ A^. 



Now we solve the system of ordinary differential equations {(C.8), (C.15)-(C.18)} numer- 
ically. For the high energy scale we take the extreme case Ai ~ 10^^ GeV and we flow to 
the weak scale, m ^ 10^ GeV, in which case t = ln{fi/A4) takes values approximately in the 
range > t > —30. 



In the SM the strong coupling (in the MS renormalization scheme for concreteness) at the Z 
pole mass is known to be PU] gs{Mz) = \J ATTas{Mz) ' ~ 1.22 and runs to ~ 0.5 at the scale 
of grand unification in a typical SU{b) model with a desert between the two scales. Thus for 
illustrative purposes we take the initial condition g{M.) = 0.5. 

The Yukawa coupling in our toy model should be thought of as the case in which the top 
quark couples to the "Higgs-like" scalar 0i. In the model studied in the main text, the ^44- 
triplet Higgs field $ is assumed to give mass to both the charged leptons and the quarks, so 
the situation y ~ ~ 1 is the one appropriate for the present study. Thus we take y{M.) = 1. 



For the scalar interactions, the purpose is to study the case Ai(A^) ~ \2{.M) ~ 1 an^ 
|Ai2(A^)| < 1. To be concrete we take Ai(A^) = A2(-M) = 1 and X^iM) = lO^^. 



The running couplings 5f(^), Ai(/i), A2(/i), and Ai2(/i) are displayed in Figs. |l(a)j |l(b) 



2(a) |2(b) , andkll respectively. The gauge coupling blows up in the deep infrared, as expected. 



and carries the Yukawa coupling with it into the nonperturbative regime. Both couplings 
remain ~ 1 down to t ~ —25, or ij,/Ai = e* ~ 10~^^. 



The quartic self-coupling Ai remains essentially fixed to its high-scale value until about the 
same scale at which g and y become non-perturbative, and the quartic-self coupling A2 runs 
down slowly from X2{Ai) = 1 to A2(e~^^A^) ~ 0.6. The difference in running between the 
two is of course due to the different roles of 0i and 02 in the model: the gauge coupling feeds 
into the running for the Yukawa coupling, which contributes only to interactions involving 
(f)i and not 02. 

The main figure of interest is Fig. [3} which shows that the off-diagonal coupling A12 runs 
down from its high-scale value Ai2(A^) = 10~^ to ever smaller values, down to Ai2(e~^^A^) ~ 
2 X 10~^. Extrapolating these qualitative results to the model in the paper, we expect that the 

^^In a two- Higgs doublet model with Higgs doublets (pi and (/)2, the VEVs must align in the charge- 
conserving vacuum so that the photon does not obtain a mass. This requires the coefficient of the interaction 
|</)|02p to be negative. Furthermore, if lm((0i) ((/)2)) = then the coefficient of (^|(/)2)^ + h.c. should also be 
negative. [5T] In the present case, the coupling A12 stands roughly for these two interactions as well as for the 
coefficient of |0ip|^2p- The running of A12 will be slow enough such that its sign does not flip, and so if we 
extrapolate to the SU{2) x C7(l)-invariant case we can just pick the appropriate sign at the high energy scale 
M, which is taken much less than the typical GUT scale, such that electromagnetism remains unbroken at 
low energy. 
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rather mild separation of scales between the weak scale and the "heavy" mass scales of either 
the Higgs doublet (f) [Eqs. (|lVg, (|IV^] or the Higgs triplet A [Eqs. ( |VL2| , ( |VL5| )] 

does not spoil the small initial values assumed for the quartic self-couplings in the scalar 
potential, which would mix the "charged lepton" and "neutrino" sectors of the model and 
introduce the sequestering problem. 
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